Single-Site Entanglement of Fermions at a Quantum Phase Transition 
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We show that the single-site entanglement of a generic spin- 1/2 fermionic lattice system can be 
used as a reliable marker of a finite-order quantum phase transition, given certain provisos. We 
discuss the information contained in the single-site entanglement measure, and provide illustrations 
from the Mott-Ifubbard metal- insulator transitions of the one- dimensional (ID) Hubbard model, 
and the (ID) Hubbard model with long-range hopping. 
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Introduction. The study of entanglement properties 
of many-particle systems has become a subject of intense 
interest. Much of the motivation comes from quantum 
information theory where entanglement is made the key 
physical resource for a variety of information process- 
ing tasks 1]. In recent work it has been suggested that 
this resource may be efficiently extracted from a solid, 
or from some other many-particle system, by scattering 
particles off the system Thermodynamic properties 
of solids have also been shown to be crucially influenced 
by entanglement properties of their microscopic degrees 
of freedom Moreover, a rapidly growing body of 
results 0, 1^ suggests that a properly chosen measure 
of entanglement may serve as a precise and convenient 
marker of a (zero-temperature) quantum phase transition 
(QPT) in a many-particle system For spin- 1/2 sys- 

tems (lattices of localized coupled qubits) a discontinuity 
(divergence) in the (derivative of the) ground state con- 
currence has been shown to be associated with a first 
(second) order QPT P| (where concurrence 7] measures 
the entanglement of two qubits selected at neighboring 
sites). For itinerant particles the picture is less clear, 
as the results here appear to depend on the choice of 
model or on the perturbation driving the transition. A 
case in point is the single-site entanglement of the one- 
dimensional (ID) Hubbard model. This measure, which 
is given by the von Neumann entropy at a single lattice 
site Q, reaches a maximum at a metal-insulator transi- 
tion driven by a change of the on-site interaction ^ . In 
contrast, the single-site entanglement diverges when one 
drives the transition by tuning the chemical potential . 

One should here realize that an onset of non-analyticity 
in a local entanglement measure is indeed expected at 
a QPT. By definition, a QPT is a point of non-analyticity 
in the ground state energy of a quantum system (caused 
by a level crossing, or, an avoided level crossing in the 
thermodynamic limit) [ll] |. Given that the elements of 
the reduced density matrix — upon which any local en- 
tanglement measure is built — are linked to the ground 
state energy, the defining non-analyticity of a QPT will 
infect also the local entanglement measure (of which 
single-site entanglement [ij, concurrence 't], and nega- 



tivity are some of the most commonly used). The 
recent proof that any entanglement measure can be ex- 
panded as a unique functional of the first derivatives of 
the ground state energy (with respect to the parameters 
that control the QPT) puts this intuition on firm ground 

The connection between entanglement and QPTs can 
also be cast in the language of statistical mechanics, 
as pointed out recently by Campos Venuti et al. 
As an example, consider the Hamiltonian density 7i(g) 
of a system that undergoes a continuous second-order 
QPT when changing a parameter g: Ti{g) — Tio -\- gA. 
Differentiating the energy density {4'o\'H{g)\ijjo) of the 
ground state | V'o) with respect to g, its singular part 
Og ^ [{ipo\A\ipo) — regular terms] will behave as Og ~ 
sgn(5 — gc)\g — gc^ as g approaches gc, implying a diver- 
gence oidOg/dg ^ \g — gc\''~^ at criticality. The singular 
term Og enters every reduced density matrix that con- 
tains a site where the operator A is defined, and it follows 
that any entanglement measure constructed from such a 
density matrix exhibits a singularity with an exponent 
related to p (barring accidental cancellations). 

Having established this linkage, one may ask how it 
can be exploited for a specific problem. For example, in 
the case of a continuous second (or higher) order QPT, 
is it possible to "read off" the critical exponent p from 
the singularity of the entanglement measure? Conversely, 
is the information provided by the singular behavior of 
a local entanglement measure already contained in the 
scaling of observables — as predicted within the usual 
statistical mechanics framework? 

In this article we address these questions by studying 
the single-site entanglement of a generic fermionic lat- 
tice system. We do so by constructing and analyzing its 
explicit representation using the Hellman-Feynman the- 
orem. We find that the single-site entanglement mea- 
sure can be used as reliable marker of a finite-order QPT 
(given certain provisos) and that it contains unique and 
useful information about the transition. The questions 
raised above will both turn out to have negative an- 
swers. As illustrations we use our construction to obtain 
the single-site entanglement at the Mott-Hubbard metal- 
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insulator transitions of the ID Hubbard model [1 
the ID Hubbard model with long-range hopping |16| . ex- 
ploiting exact results for the ground state properties of 
these models. We stress that our analysis can be easily 
adapted so as to apply to a system of localized spins, 
with no change in the general results. Specifically, the 
questions raised above are answered in the negative also 
for coupled qubit (spin-1/2) systems. Our reason for fo- 
cusing on fermionic systems is simply that these are less 
well understood. With our contribution we hope to dis- 
pel some of the perceived difhculties attached to their 
treatment. 

Single-site entanglement and QPTs. Let us first 
recall that the concept of quantum entanglement of in- 
distinguishable fermions [bosons] suffers from a certain 
ambiguity since the accessible state space contains only 
antisymmetrized [symmetrized] states and hence lacks a 
direct product structure. The simplest way around this 
problem is to use an occupation number representation 
UJ. For spin-1/2 fermions one thus takes | n)j =| Q)j, 
IT)i) and|ti)j as local basis states, with j = 1,2, L 
indexing the corresponding lattice sites. In this way the 
product structure of the state space is manifestly recov- 
ered, with the representation spanned by the 4^ basis 
states I n)i(g) | 71)2 (8) ...(§) | n)^. One may now proceed 
as usual and partition the system into two parts A and 
B, with the entanglement (yon Neumann) entropy f of a 
pure state | ip) defined by 0] 

£ = -TripAlog^PA). (1) 

The reduced density matrix pA is calculated from the 
full density matrix p =\ip){'ip \hy taking the trace over 
the local states belonging to B: pA = Ttb{p)- By choos- 
ing A as a single site (assuming translational invariance) 
with B the rest of the system, one obtains the single-site 
entanglement. One should note that in the occupation 
number representation the subsystems A and B corre- 
spond to fermionic modes (empty sites, singly occupied 
sites with spin up or down, doubly occupied sites) and 
not to particles. In this sense the notion of fermionic 
(and similarly, bosonic) entanglement is different from 
the text book example with spatially separated particles. 

Given the occupation number representation it is 
straightforward to verify that the reduced ground state 
density matrix pj for a single site j is diagonal, provided 
that the ground state \ ipo) is a superposition of basis 
states with the same number of particles and the same 
total spin. Introducing the ground state expectation val- 
ues for a single site to be doubly occupied (102), singly 
occupied by a fermion with spin-up [spin-down], (u;|[j]), 
or empty (wq), and assuming that the system is transla- 
tionally invariant, we write: 



and where in Eq. (0) 



cj^Cjo- is the number operator 



W2 
Wq 



{■)po\hj^hji\^o) 
(V'ol^iTlV'o) - u 

(V'ol^-ilV^o) - W2 



{ipo\nj^\ipa} -~W2 = -+m-W2 
n 



m — W2 



that samples site j for a fermion of spin a =T,i, n = 
('0o|fT-j| + fij^l^po) is the average single site occupation in 
the ground state, and m = {l/2){tpQ\nj'^ — nji\ipo) is the 
ground state magnetization per site. It follows that 



Pj 



E 

a=0,T,i 



Wa\a)j{a\j +W2\U)j{U\: 



(3) 



Combining Eqs. |5J, and the single-site entan- 
glement takes the form 



W2 log2 



W2 



!l-m-W2)\og2{^ 



m — W2 

(1 - n -I- W2) log2 (1 - n -I- W2) ■ 



W2 l0g2 W2 

(4) 



Let us now consider a fermion system with Hamiltonian 
density Ti-ig) = TLo + gh. that exhibits a QPT for 
some value gc of g (with A the conjugate operator, 
and with all other control parameters kept fixed and 
absorbed as part of T-Lq). By definition, a QPT of 
^th Qj-jjgj. implies a divergence or a discontinuity in 
the k^^ derivative d^eo/dg^ of the ground state en- 
ergy density eo = ('0o|'W(5)|V'o), with all derivatives 
of order < k being finite and continuous. Defining 
Og = [(?/'o|A|?/'o) — regular terms] (equal to [deo/dg 
- regular terms] by the Hellman-Feynman theorem), 
it follows that d*'~^Og/dg*'~^ has a divergence or a 
discontinuity at g = g,.. With these preliminaries we can 
now prove the following 



Proposition 

Consider a spin-1/2 translationally invariant fermionic 
system with a Hamiltonian density Ti.{g) = Tio + gA 
that conserves particle number and total spin, and where 
Og = [{ipo\M''Po) ~ regular terms] is a linear combination 
of m,n and/or 'W2. It follows that a divergence or a dis- 
continuity in the (k — 1)^^ derivative of the single-site 
entanglement with respect to g (with all derivatives of 
order < fc—1 being finite and continuous) signals that the 
system undergoes a k^^ order QPT. 

Proof 

The proof is elementary. Repeated differentiation of Eq. 
(O yields 



I — n + W2 



(2) 
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n 



Q J— j- log2 (^2) + — j log2 (1 — n + W2) + terms containing lower-order derivatives. (5) 



By assumption all derivatives with respect to g of order 
< k —\ are finite and continuous. Any singularity in 
Qk-igjQgk-i j^^gi; hence reside in terms containing 
derivatives of order fc — 1 . Since Og is a linear combina- 
tion of m, n and W2, the proposition follows. I 

Several comments are in order. First note that the 
constraint that Og should be some linear combination of 
m, n and/or W2 is much less restrictive than may first ap- 
pear to be the case. In fact, for a generic fermionic QPT 
caused by a change of an interaction or an external per- 
turbation that couples only to single sites, Og is identical 
to W2 (with the transition driven by an on-site fermion- 
fermion interaction, g = u), m (with the transition driven 
by a magnetic field, g = h), or n (with the transition 
driven by a chemical potential, g = One may think 
that the tight link between the scaling of d^^^E /dg^^^ 
and that of d^^^Og/dg^^^ would allow for the critical 
exponent that controls Og to be immediately extracted 
from d''~^£/dg''~^. This is not so, however. As an ex- 
ample, take a second order QPT (fc — 2) with Og ~ W2, 
where dw2/du ^ \u ~ Uc\'^~^ — )■ cxd as g ^ gc = Uc- By 
inspection of Eq. © one then notes that the leading 
scaling of d£ / dg will be governed by the same exponent 
p only if m and n are independent of W2, or, depend 
on W2 as a power with exponent > 1. Whether this is 
the case typically requires that one has access to an ex- 
act solution of the model, and in any event can only be 
determined on a case-to-case basis. Turning to the log- 
arithmic factors in |SJ one realizes that these will cause 
logarithmic divergences if one or several of the occupa- 
tion parameters •wo,w^,wi,W2 vanish at the transition 
(cf. the parameterization in ij^J). Such logarithmic cor- 
rections, multiplying the leading scaling of d'^~^£ /dg'^^^ 
inherited from Og, thus signal a change of the dimen- 
sion of the accessible local Hilbert space as the system 
undergoes the transition. This is a useful and impor- 
tant property of the single-site entanglement scaling not 
shared by the scaling of Og or its derivatives. One should 
here note that a spurious signaling of a fc^'^ order QPT 
by a divergence in d''~^£ /dg'''~^ caused by a vanishing 
occupation parameter is blocked by the constraint in the 
proposition that all lower-order derivatives of£ are finite. 
(Although maybe hard to realize, one may envision a sys- 
tem where one or several local basis states get excluded 
when tuning some parameter in the Hamiltonian [imply- 
ing the vanishing of an occupation parameter] without 



the occurrence of a QPT.) 

Using the diagnostics supplied by our proposition, are 
we guaranteed to catch all fermionic QPTs? The an- 
swer is negative. First, the diagnostics obviously fails 
for a QPT of infinite order [l^ . a Berezinskii-Kosterlitz- 
Thouless (BKT) type transition being a case in point 
|l9l|. Secondly, and more insidious, a system may exhibit 
a QPT of finite order, but with the single-site entangle- 
ment and its derivatives still remaining regular. This 
happens if all local basis states \ n)j =|0)j, |t)j, and 
Ti)i become equally populated as one approaches the 
transition. As seen from 0, the (fc — 1)^^ derivative 
terms then vanish identically, killing the signal of the 
QPT. The simultaneous vanishing of d£/dg implies that 
£ has a local extremum at the transition (expected to be 
a maximum since in this case all local basis states are 
equally represented in the make-up of the many-particle 
ground state). However, one cannot a priori exclude that 
£ is at an extremum without the occurrence of a QPT. 
Hence, an extremum of the single-site entanglement does 
not necessarily signal a QPT. Whether a QPT is present 
or not in this case requires information beyond that pro- 
vided by the entanglement measure. 

Having exposed the general features of entanglement 
scaling at a fermionic QPT, let us look at two examples. 

Case studies. Consider first the ordinary ID Hub- 
bard model 

L L 

H ^ - ^ (cj^Cj+io- h.c.) fii^fiii (6) 

i=l 1=1 

with the first term describing hopping of electrons be- 
tween neighboring sites, and with the second term an 
effective on-site interaction of strength u. At half-filling 
of the lattice, n = 1, the model exhibits a QPT at u = 0, 
separating a Mott insulating phase (u > 0) from a metal- 
lic phase {u < 0). The ground state energy density be- 
comes non-analytic at the transition, but allows for an 
asymptotic power series expansion with all derivatives 
being finite and continuous jl^]. The QPT is thus of 
infinite order, and can be shown to belong to the BKT 
universality class As found by Gu et al, the single- 
site entanglement has a maximum at the transition. This 
reflects the equipartition of empty-, singly- and dou- 
bly occupied local states when u — (non-interacting 
fermions). The transition is thus special on two counts: 
it is of infinite order and it supports an equipartition of 
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local states. This makes it an exceptional example of 
a fermionic QPT, where no information can be deduced 
from the entanglement measure. 

A metal-insulator transition can also be triggered when 
M > by connecting the system to a particle reservoir and 
tuning the chemical potential g = fi: When n < 1 the 
system is metallic, but turns into an insulator at the crit- 
ical pointiic = 2 — 4 Ji(w)[a;(l-|-exp(cj'u/2))]~^ where 
n = 1 The transition is second order with a diver- 
gent charge susceptibility Xc = dn/d^ ~ — /Xc|~^^^. As 
shown in Ref. 0, the derivative of the critical single-site 
entanglement for finite u is precisely given by Xc up to a 
multiplicative constant: = —C{u)xc- In the limit 

M — > oo the empty local states get suppressed at the tran- 
sition and the scaling of dS/dfx picks up a logarithmic 
correction 0]: d£/dpL = Xc(ln|^ — /^d + const. )/21n2. 
Both behaviors well illustrate our general discussion 
above: For finite u the logarithms in Eq. jSJ add up 
to the M— dependent constant C(m), whereas in the limit 
M — > oo the entanglement measure detects a change in 
the dimension of the local Hilbert space, as signaled by 
the logarithmic correction to the leading scaling. 

As a second example, let us consider the ID Hubbard 
model with long- rang e hopping, introduced by Gebhard 
and Ruckenstein |1(tI |: 

L L 
= ^ tlmc\^Cma ^ u'^h^'ihtl (7) 

<T=Ta 

with = i(— — to)^-^. The ground state en- 
ergy density at half- filling is given by eo — {un — Uc(l — 
n)n)/4- (1/(24mMc))[(u + Mc)^ - iiu + u^f - 'iuucnf'^] 
with Uc = 27r the critical point This implies that 

W2 = deo/du has a discontinuity in its second order 
derivative with respect to u at Uc and hence the tran- 
sition is third order. From Eq. Q with n = 1 it fol- 
lows that the single site entanglement can be written as 
£ = —(1 — 2w2)log2(l/2 — W2) — 2w2log2(w2) when no 
magnetic field is present (i.e. to = 0), and one imme- 
diately verifies that d^£/du^ is also discontinuous at the 
transition point Uc- Since the local basis states do not be- 
come equally populated at Uc — in contrast to the u = 
metal-insulator transition of the ordinary Hubbard model 
— the single-site entanglement here provides an accurate 
diagnostics of the transition. 

One can also drive a Mott-Hubbard metal-insulator 
transition by tuning the chemical potential when u > Uc, 
in exact analogy with the ordinary Hubbard model. Ex- 
pressing n as a function of /z, and applying the Hellman- 
Feynman theorem to the ground state energy eo above, 
one obtains a discontinuity in dn/d^ at /i = /ic = |22l |. 
Eq. (jsj immediately implies that d£ /d^ is also discontin- 
uous at fj, — lie, with the transition being second order. 
In the limit u 00 this discontinuity is multiplied by 
a logarithmic divergent factor when ^ — > /Zc-, reflecting 
the suppression of empty states in this case. 



Summary. We have shown that a generic finite-order 
quantum phase transition in a spin-1/2 fermionic lattice 
system can be consistently identified and characterized 
by studying the behavior of the single-site entanglement 
and its derivatives with respect to the parameter that 
controls the transition. Extensions to cases where the 
transition is driven by an interaction or a field that cou- 
ples to pairs or clusters of lattice sites (like the extended 
Hubbard model 0|) is conceptually straightforward, al- 
beit technically more demanding. We hope to return to 
this problem in a future publication. 
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